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PACS 74.45.+C - Proximity effects; Andreev reflection; SN and SNS junctions 
PACS 73.23.Hk - Coulomb blockade; single-electron tunneling 
PACS 72 . 25 . -b - Spin polarized transport 

Abstract. - We investigate electronic transport in a three-terminal hybrid system, composed 
by an interacting quantum dot tunnel coupled to one superconducting, one ferromagnetic, and 
one normal lead. Despite the tendency of the charging energy to suppress the superconducting 
proximity effect when the quantum dot is in equilibrium, the non-equilibrium proximity effect can 
give rise to a large Andreev current. The presence of the ferromagnet can lead to a finite spin 
accumulation on the dot. We find that the interplay of the Andreev current and spin accumulation 
can generate a pure spin current, with no associated charge transport, in the normal lead. This 
situation is realised by tuning the quantum-dot spectrum by means of a gate voltage. 
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Introduction. — Sub-gap transport through hybrid 
superconducting-normal interfaces is sustained by An- 
dreev reflection [1,2 . In an Andreev-reflection process, 
two electrons vi^ith opposite spins enter the superconduc- 
tor from the normal region forming a Cooper pair. In 
the reverse process, a Cooper pair breaks in the super- 
conductor and its constituents, two electrons with oppo- 
site spins, enter the normal region. In quantum dots An- 
dreev reflection has been studied theoretically including 
Coulomb interaction and non-equilibrium [5H5]. Andreev 
reflection in quantum dots couples the electron and hole 
dynamics and it leads to the formation of Andreev bound 
states. In the limit of infinite superconducting gap, when 
quasi-particle transport is inaccessible, the coupling be- 
tween quantum dot and superconductor can be treated 
exactly [9HT4). Recently, several experiments have mea- 
sured electronic transport through quantum dots coupled 
to superconducting leads [TSHSn]- 

Hybrid systems containing ferromagnetic materials are 
of interest for the possibilities they open up to gener- 
ate and control spin currents. The current between two 
tunnel-coupled ferromagnets depends strongly on the rel- 
ative magnetization direction of the two ferromagnets. 
For parallel alignments the tunneling current is enhanced, 
while for antiparallel alignments the current is suppressed. 
This effect is known as the tunneling magneto resistance 



(TMR) [5d;. 

In quantum dots connected to ferromagnets, spin accumu- 
lation in the dot plays an important role. Theoretically, 
for such systems complex transport properties as negative 
differential conductance [311132] and spin precession [33134] 
have been predicted. Experimentally, quantum dots real- 
ized in carbon nanotubes have already been coupled to 
ferromagnets [551 - 157] . 

Cooper pairs in a BCS superconductor are made up of a 
spin up and spin down electron in a singlet state. Hence 
the combination of superconducting and ferromagnetic 
materials in a nano-structure is expected to give rise to 
rich spin physics and to non-local effects due to the en- 
tanglement intrinsic to the Cooper pairs. A good example 
is crossed Andreev reflection in three terminal setups con- 
sisting of a superconductor and two ferromagnets [38H42] . 

Recently, S. Das et al. have proposed a scheme for pure 
spin currents in a three-terminal setup, where a ferromag- 
net, a superconductor, and a normal conductor are con- 
nected via quantum wires [43]. In the present work, we 
consider a similar setup but involving a quantum dot. In 
particular, we consider an interacting, single-level quan- 
tum dot tunnel coupled to a superconducting, a ferromag- 
netic, and a normal lead, see Fig. [T] We exploit the tun- 
ability of the quantum-dot spectrum and non-equilibrium 
to generate a pure spin current in the normal lead, i.e. a 
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Fig. 1: (color online) Setup: a quantum dot is tunnel coupled 
to a superconducting, a ferromagnetic, and a normal lead. 



spin current with no charge current. 

Model. — The Hamiltonian of the three-terminal sys- 
tem under investigation can be written as H ^ H^ot + 
Sr)(^') + -fftunn,))), whcre the different leads are identified 
by the label 77 — N,F,S corresponding to normal, fer- 
romagnetic and superconducting, respectively. The leads 
are described by the Hamiltonians 

-^"77 = ^ £rjkacl^.„Cr,kcr - Sjj.S ^ (Ac^_fe4_C^fe-|- -|- H.C.) , 
ka k 

(1) 

where the single-particle energies Cr^ka are spin dependent 
only for the ferromagnet, and the BCS pair-potential A is 
nonvanishing only for the superconductor. The annihila- 
tion and creation operators for electrons in lead 77 are c^^ka 
and c^j,^, respectively. Superconductivity is described by 
the mean-field BCS Hamiltonian with a pair potential A, 
which has been chosen to be real. The chemical potential 
of the superconductor is taken as reference for the ener- 
gies, i.e. ^s = 0. 

The single-level quantum dot is described by the 
Anderson-impurity model: 



Hdot = ^ <id)ada + Un^ni, 



where e is the energy of the spin-degenerate single-particle 
level, and U the on-site Coulomb repulsion. The opera- 
tors dj annihilate (create) a dot electron with spin a, and 
Tier = dl-da is the corresponding number operator. Elec- 
tron tunneling between dot and leads is modeled by the 
tunneling Hamiltonians 



H< 



tunn,?7 



y.T. 



Sfeo- " 



H.C. 



(3) 



where for simplicity the tunnel-matrix elements Vri are 
assumed to be independent both of spin and momen- 
tum. The strength of the tunnel coupling is captured 
by the spin-resolved tunnel-coupling strengths defined as 
Trier = 27r|K,p J2k ^i^ - ^vkcr) = 27r|l/^ ^p^^r , with pna- be- 
ing the density of states of electrons with spin a in lead 
77, which we assume to be constant in the energy win- 
dow relevant for transport. The mean level broadening 



is then defined as F,, — ^^^r^vc- "^^^ spin asymme- 
try of the density of states in the ferromagnet at the 
Fermi energy is characterized by the spin polarization 
P = {PFt - PFi)/{pF\ + PFl)- 

Method. — To obtain an effective description of the 
system, we integrate out the leads' degrees of freedom. 
For voltages and temperatures smaller than the supercon- 
ducting gap A, quasiparticle excitations in the supercon- 
ductor are suppressed. In the following, we consider the 
limit A — > 00, for which they are completely inaccessible 
and only subgap-transport is possible from or to the su- 
perconductor. This limit implies that A also exceeds the 
bias voltages applied to the leads as well as the charging 
energy U of the dot (such that the energies for charge ex- 
citations in the dot lie within the superconducting energy 
gap). The latter condition is experimentally challenging 
but not impossible, see e.g. Ref. [TB]. In the limit A -^ 00, 
one can integrate out the coupling to the superconductor, 
which leads to an effective Hamiltonian [9l-[12] describing 
the proximizcd quantum dot: 



^cff — 2_^ ^' 



d^da + Un-^ni Trdldl Trd^d- 



(4) 



The Hilbert space is spanned by the four eigenstates \x) G 
(h) Jt) , li> , 1+)}, with eigenenergies E_ = ^- €a, E-f 



n 



E, 



e, and E-l 



EA, where |t) and \i) correspond 



to a singly occupied dot and 



1 

1 

V2 



^ + ^'«> 



i'«^ 



2eA 



\d) 



^^^'^^ 



are superpositions of an empty, |0), and doubly occupied 
(2) dot, \d) = rfldj |0) . Here, S = 2e+U is the detuning and 



4 



£A is given by e^ = Y (e + t)^ + 
Due to tunnel coupling to a superconductor the particle 
and hole sector of the Hilbert space are mixed, leading to 
four Andreev bound-state energies: 



EA,Y,f = 7 -^+1£A, 



(5) 



with 7', 7 = ±. The Andreev bound-state energies 
are the excitation energies of the proximized dot for 
transitions between the singly occupied states |t) or \\.) 
and the superposition states |— ) or |+). We consider 
the regime \S\ < \/U'^ — F^, for which the inequality 
Ea-- < Ea-,+ < < Ea,+ - < Ea,+,+ holds. 



The dot is described by its reduced density matrix 
Prcd, whose matrix elements are P^^^ = (xi|yOred|X2)- 
We consider only the case of weak coupling to the non- 
superconducting leads, F7v,Fi? < fc^T and we compute 
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results up to first order in these coupfings. We also as- 
sume that FatjEf <C tA- When these conditions are ful- 
filled, we need to retain only the diagonal elements of the 
density matrix in the basis of the eigenstates of eq. Q. 
The diagonal elements are the occupation probabilities of 
the corresponding states and we denote them by P^ = P^. 
In the stationary limit, the occupation probabilities 
obey the master equation 



E 



iW^x'Px' 



WrxPx) 



0, 



(6) 



where W-^n-^i is the transition rate from state x' to state 
x", which can be computed by means of Fermi's golden 
rule. The current of electrons with spin a out of lead 77 
can be expressed as 



x^x' 



W''"" P 



(7) 



where W27,, = Vl^ sW!'^, and VF''^ is the sum of all rates 

X X — ^ XX XX 

that describe transitions from x to x in which in total s 
electrons of spin a are removed from lead rj. The charge 
current and the spin current out of lead rj are given, re- 
spectively, by 






tS - 



(8) 
(9) 



where cq = e and cs = hj^. 



Results. — We discuss the case of equal tunnel- 
coupling strengths to the normal and ferromagnetic leads, 
Ftv — Fi? = F < fcsT , and focus on the regime of positive 
bias, /^AT > 0. The case [itq < is obtained from the sym- 
metry transformation /ijv — >■ —p-N, \J^f — > — Mf, <5 — >■ —J, 
and Jf^ — > —Jjq- 

The spin current in the normal lead 



J j^ — — OzF^v 



{fU^^A 



(10) 



:)^!^{Ea^-A\ (1-^ 



is proportional to the spin accumulation Sz — 
(Ft - n) /2, where f+{Lu) = /jv(c^) = [1 + cxp{^^)]-^ 
is the Fermi function of the normal lead and f^{uj) — 
1 — /Ar(a;). Ineq. ([TU]) , the effect of the voltage bias applied 
to the ferromagnetic lead, fip, is contained in the spin ac- 
cumulation Sz- For small chemical potentials applied to 
the normal lead, Ea,-,+ < ^J,N < Ea,+.-, the spin current 
is exponentially suppressed. For intermediate chemical 
potentials of the normal lead, Ea,+.- < IJ-n < Ea.+,+, 
and arbitrary hf, the spin current simplifies to 






'SzT 1 



S 
2^ 



(11) 



In the same bias regime, the charge current is given by the 
expression 



j: 



N 



F 



1- 



5 

2^A 



Pi 



1 



S 
2^ 



P- - 



-P+ 



^A 



(12) 

where Pi = P-^ + P^ is the probability for the dot to be 
singly occupied. We note that the charge current in inter- 
mediate bias regime can be tuned to zero by means of the 
dot-level position e. 

On the other hand, in the large-bias regime for the normal 
lead, Ea,+,+ < fJ-N, where processes with electrons enter- 
ing the normal conductor are exponentially suppressed, 
there will always be a finite Andreev charge current flow- 
ing out of the normal lead into the superconductor, i.e., 
the charge current will be finite for any value of the level 
position. 

In conclusion, a pure spin current is only possible in the 
intermediate-bias regime, which we will focus on in the 
following. 

How to tune the system parameters to obtain a pure 
spin current is illustrated in Fig. [21 where the charge and 
spin currents are plotted as a function of the level position 
(tunable by the gate voltage) and of the chemical potential 
of the ferromagnetic lead ^f (tunable by the bias voltage). 
For the value of the level position corresponding to the thin 
line, the charge current vanishes for Ea,+,+ < ^^f, while 
the spin current remains finite. The level position eg that 
tunes the charge current in the normal lead to zero, is a 
function of the charging energy U , the polarization p, and 
the tunnel-coupling strength to the superconductor F5: 




-p4 + 27p2 - 18 + (9 - p2) ^pi - 36p2 + 36 



l-p2 



(13) 



The equation is valid for the normal lead being in 
the intermediate-bias regime, Ea,-,-, Ea,-,+ , Ea,+.- < 
fiN < Ea.+,+ , and the ferromagnetic lead being in the 
large-bias regime, Ea,+,+ < Hf- For its derivation we have 
assumed that all occurring Fermi functions can be approx- 
imated to either zero or one, i.e. that l^N — EA.-y'.-yl, \^J'F — 

Ea,j',j\ ^ ksT for all possible values of 7 and 7'. 

Increasing the polarization p increases the spin accumu- 
lation on the dot and reduces the proximization of the dot 
by the superconductor. Thus the probability for the dot 
to be singly occupied. Pi, is increased, while the proba- 
bility for the dot to be in state |+), P+, is reduced. In 
order to keep the charge current in the normal lead tuned 
to zero for increased polarizations the prefactor of Pi , i.e. 
(1 — (5/2e^), in eq. (fT2|) has to be reduced by increasing 
e = Eq. 

For moderate values of the polarization p the monotoni- 
cally increasing function £o{p) grows slowly. Only for large 
polarizations ea{p) grows fast. 
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Fig. 2: (Color online) Density plot of a) the charge current and 
b) the spin current in the normal lead as a function of the level 
position e and the chemical potential of the ferromagnetic lead 
fjLF- The thin line indicates the value of the level position for 
which the charge current vanishes in the high-bias regime for 
the ferromagnetic lead. The other parameters are /ijv = 0.5L*', 
Ts = 0.2(7, r = O.OOlt/, p = 0.6, and keT = 0.0117. 



When eo becomes of the same order of the single-particle 
level spacing in the quantum dot, the single-level model 
for the dot ceases to be valid. However, for the system 
under consideration this situation is not realized for real- 
istic polarizations of the ferromagnetic lead. In fact, for 
Ts = 0.2U and p = 0.995 we find that the level position 
that tunes the charge current in the normal lead to zero 
only reaches eo ~ 0. 

The pure spin current for e = eo is independent of 
the tunnel-coupling strength to the superconductor Ts 
(though a finite Ts is required to establish the spin cur- 
rent and the Andreev current tuning the charge current 
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Fig. 3: (color online) Spin current J^ and charge current J^ 
as a function of the chemical potential of the ferromagnetic 
lead ^iF for e = e,,, Miv = 0.517, Fs = 0.217, F = 0.00117, p = 0.6, 
and ksT = O.OIC/. 



to zero). 

Figure [3] shows the spin and charge currents in the nor- 
mal lead for e = eo, as a function of the chemical po- 
tential of the ferromagnetic lead fip- The dependence of 
the pure spin current on the polarization of the ferromag- 
net is shown in Fig. 21 With increasing polarization the 
spin current increases until a maximum is reached. For 
stronger polarizations, the charge current flowing out of 
the ferromagnet is reduced due to the spin accumulation, 
leading to reduction of the spin current in the normal lead. 

The occurrence of a pure spin current can be under- 
stood examining the transport processes between dot and 
ferromagnetic / normal lead. Since the ferromagnetic and 
normal leads are only weakly coupled to the quantum dot, 
corresponding currents are carried by sequential electron 
tunneling. 

The possible transport processes between dot and ferro- 
magnetic lead in the large bias regime, i.e. Ea,+,-{- < /^f, 
are illustrated in Figure [S] a). The processes involving mi- 
nority spins, depicted by shorter arrows, are suppressed. 
In all the processes electrons are transferred from the fer- 
romagnetic lead to the dot. All four dot states can be 
occupied. Due to the spin polarization of the ferromag- 
netic lead, an imbalance in the occupation of the dot states 
It) and It), and thus a finite spin accumulation Sz, is gen- 
erated. 

The possible transport processes between dot and normal 
lead in the intermediate bias regime, i.e. EA,-i-,- < jJ-n < 
Ea.+,-^i Eire shown in Figure [5]b). Because of the spin ac- 
cumulation on the dot, i.e. P^ > Pi, transitions from state 
It) to state |— ) take place more frequently than transitions 
from It) to |— ). This leads to a finite spin current in the 
normal lead. Black arrows represent processes where an 
electron flows out of the normal lead, while white arrows 
indicate processes where electrons enter the normal lead. 
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Fig. 4: (color online) Spin current in the normal lead as a func- 
tion of the polarization of the ferromagnet p for e = eo, Miv = 
0.5i7, Vs = 0.2U, r = O.OOlf/, and kgT = O.OIU. 
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Fig. 5: Possible tunneling processes between dot and a) ferro- 
magnetic lead with large chemical potential Ea,+,+ < fJ,F and 
b) normal lead with intermediate chemical potential Ea,+,- < 
j-iN < Ea,+,+'- the arrow direction depicts transitions of dot 
states. With black (white) arrows an electron leaves (enters) 
the corresponding lead. Tunneling of minority spins (short ar- 
rows) is suppressed. 



The latter processes require the dot to be in state |+). 
The occupation of dot state |-|-) is provided by the fer- 
romagnet. By controlling the gate voltage, and thus the 
occupation of the different dot states, the charge current 
in the normal lead can be tuned to zero yielding a pure 
spin current. This pure spin current is generated by the 
interplay of superconductivity, non-equilibrium, and spin- 
dependent transport: the superconducting lead proximizes 
the quantum dot, leading to the superposition states |— ) 
and |-|-). The biased ferromagnetic lead first provides oc- 
cupation of the dot state |-|-), which gives the possibility 
to tune the charge current in the normal to zero. Second, 
it causes spins to accumulate on the dot inducing a spin 
current in the normal lead. 

Conclusions. — Spin-dependent transport in a three- 
teminal structure consisting of a quantum dot tunnel cou- 
pled to a superconducting, a ferromagnetic and a normal 
lead has been investigated in the large-superconducting- 
gap regime, by means of a master-equation approach. It 
has been established that such a system can be used to 
generate a pure spin current in the normal lead, exploit- 
ing the tunability of the dot's spectrum, non-equilibrium 
driven by finite bias voltages, spin accumulation due to 
the ferromagnet, and the Andreev bound states induced 
by the superconductor. 
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